:arXiv:cond-mat/9506101v2 23 Jun 1995; 



o 

cj 
M 



TABLES 



2 



WIS-95-ED 

Lack of Self- Averaging in Critical Disordered Systems 



Shai Wiseman and Eytan Domany 
Department of Physics of Complex Systems, Weizmann Institute of science, 

Rehovot 76100 Israel 
(June 18, 1995) 



Abstract 



We consider the sample to sample fluctuations that occur in the value of a 
thermodynamic quantity P in an ensemble of finite systems with quenched 
disorder, at equilibrium. The variance of P, Vp, which characterizes these 
fluctuations is calculated as a function of the systems' linear size I, focusing 
on the behavior at the critical point. The specific model considered is the 
bond-disordered Ashkin- Teller model on a square lattice. Using Monte Carlo 
simulations, several bond-disordered Ashkin- Teller models were examined, in- 
cluding the bond-disordered Ising model and the bond-disordered four-state 
Potts model. It was found that far from criticality all thermodynamic quanti- 
ties which were examined (energy, magnetization, specific heat, susceptibility) 
are strongly self averaging, that is Vp ~ l~ d (where d = 2 is the dimension). 
At criticality though, the results indicate that the magnetization M and the 
susceptibility \ are non ser f averaging, i.e. -p-, ^ -/-> 0. The energy E at 
criticality is clearly weakly self averaging, that is Ve ~ l~ Vv with < y v < d. 
Less conclusively, and possibly only as a transient behavior, the specific heat 
too is found to be weakly self averaging. A phenomenological theory of finite 
size scaling for disordered systems is developed, based on physical consider- 
ations similar to those leading to the Harris criterion. Its main prediction 
is that when the specific heat exponent a < (a of the disordered model) 
then, for a quantity P which scales as l p at criticality, its variance Vp will 
scale asymptotically as l 2p+ ~ . The theory is not applicable in the asymptotic 
limit (/ — ► oo) to the bond-disordered Ashkin- Teller model where ^ = 0+. 
Nonetheless in the accessible range of lattice sizes we found very good agree- 
ment between the theory and the data for V x and Ve- The theory may also 
be compatible with the data for the variance of the magnetization Vm and 
the variance of the specific heat Vo, but evidence for this is less convincing. 

75.50.Lk 75.40Mg, 75.10Nr, 75.40Cx 
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I. INTRODUCTION 



How is the critical behavior affected by the introduction of disorder (usually dilution or 
bond-randomness) into a model? This question has been extensively studied Q experimen- 
tally, analytically || and numerically [[| for quite some time now. Many studies concentrate 
on finding out to which universality class certain disordered models belong, e.g. calculating 
critical exponents. In this work we consider a different aspect of the same question. The 
measurement of any density of an extensive thermodynamic property P (e.g. P = E, M, Ch 
or x) i n a disordered system may hypothetically be done in the following way. An ensemble 
of macroscopic disordered samples of size I is prepared; denote by x a sample with a partic- 
ular random realization of the quenched disorder. Now in each sample x, P x (t) is measured 
over a long time interval, and P x , the average over time t is calculated. Close to the critical 
point the measurement of P x will require long times due to large thermal fluctuations which 
will occur. In addition, since in every sample a different configuration of the quenched dis- 
order is present, a different value for P x will be measured. Next, the average of P x over the 
ensemble [P x \ ([...] stands for an ensemble average over the different samples) is calculated 
and so is its variance 



Assume that the time interval of the measurement was long enough so that thermal fluc- 
tuations in P x (t) were averaged out perfectly and P x may be considered to be exact. The 
question then rises: How will the variance Vp change as the critical temperature is ap- 
proached or as the correlation length £ [[| is increased? This question, which concerns the 
way in which disorder affects the behavior of systems near their critical point, is approached 
in this work using the framework of finite size scaling. 

A common practice in Monte Carlo (MC) simulations is to examine the critical behavior 
by simulating a system at its critical temperature T c °° and changing the lattice size I. Ac- 
cording to the theory of finite size scaling |5| the lattice size I sets the scale of the correlation 
length in such a finite system. Thus the dependence of P on £ in an infinite system close 
to criticality is substituted by dependence on I in a finite system at criticality. When a 
disordered system is considered then many samples need to be simulated in order to obtain 
estimates of P which are averaged over the disorder. In this case, the question, which is the 
main theme of this work, would be: how does the effect of disorder on the sample to sample 
fluctuations in P x change, as the lattice size I is increased at the critical temperature? Or 
how does Vp scale with Z? This question is not only of theoretical interest in its own right, 
but also of practical interest for MC studies of critical disordered systems. If the relative 
variance Vp/[P} 2 decreases with increasing I then the number of samples needed to obtain 
[P] to a given accuracy goes down with increasing I. If, on the other hand, Vp/[P] 2 is in- 
dependent of /, then the number of samples which need to be simulated is independent of / 
and the total amount of work rises very strongly with I. 

The issue which we study in this work should not be confused with two closely related 
issues. The first is usually referred to as the property of self-averaging of additive (extensive) 
quantities in disordered systems || . Consider again the ensemble of macroscopic disordered 
samples of size I. The question is then whether 



Vp = [{p x -[p x ]Y]. 



(i) 



V P /[Pf ^ as I 



oo . 



(2) 
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If so, then the measurement of P x in one very large sample x which occurs with reasonable 
probability will provide a good estimate of the ensemble average. This is very important 
for the comparison of theoretical work, where the configurational average is taken, with 
experiments, where only a large single sample is examined. As first argued by Brout J7|, we 
may divide the sample x into n large subsamples (much larger than the correlation length 
£). If we assume that the coupling between neighboring subsystems is negligible, then the 
value of any density of an extensive quantity over the whole sample is equal to the average 
of the (independent) values of this quantity over the subsamples. Provided the probability 
distribution of the P's of the subsamples has a finite variance, then according to the Central 
Limit theorem the value of P x is distributed with a Gaussian probability distribution around 
its mean [PJ. The square of the width of the Gaussian, Vp, is proportional to - ~ l~ d . In 
this case (||) is fulfilled, and P is called self-averaging. 

The quantities which are studied here are all densities of extensive self averaging quan- 
tities (far from criticality). Nonetheless, note that our question, as it was formulated for 
macroscopic samples (I ^> £), concerned the dependence of V on the correlation length £ 
and not on the sample size I. On the other hand, as we will examine finite samples of size 
I at criticality where £ ~ /, the Brout argument does not hold, since the average of P over 
neighboring subsamples may not be considered as independent. Thus at criticality there 
is no reason to expect that Vp ~ l~ d . An example for a phase transition, where sample 
to sample fluctuations result in non self averaging of certain quantities, is the percolation 
transition. It has been shown || that the resistive susceptibility and the conductivity are 
non self averaging at the percolation threshold. 

A second related issue is that of self-averaging in homogeneous systems. This question 
concerns the thermal fluctuations in the value of a density P in a homogeneous system of 
size I. Define the thermal variance as cx|> = ((P — (P)) 2 ), where (. . .) denotes thermal or 
time averaging. The following notions (slightly modified) have been introduced by Milchev 
Binder and Heermann If a T /(P) 2 — > as I — > oo then P is self averaging otherwise 



it is said to exhibit lack of self averaging. If 

4/(p) 2 ~ r d (3) 

then P is strongly self averaging. If 

a 2 T /{P) 2 ~ r Xl and < xi < d (4) 

then P is weakly self averaging. When / 3> £ it was found |||IIJ that averages of simple 
densities such as E, M are strongly self averaging while quantities obtained from the fluctua- 
tions of these densities such as the specific heat C and susceptibility x are non se h averaging. 
At criticality the singular part of the energy E is weakly self averaging while C, M and x 
exhibit lack of self averaging. For example (M) 2 ~ r 2p/u and of M ~ x/l d ~ V lv ~ d = l~ 2l3/u , 
so that M is non self averaging. 

The issues of self averaging in disordered systems and homogeneous systems concern 
the asymptotic behaviour of the fluctuations due to disorder and the thermal fluctuations 
respectively as the system size is increased. While self averaging in homogeneous systems at 
criticality has been addressed previously j this study involves the question of self aver- 
aging in disordered systems at criticality. With the increase in the available computational 
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power, a numerical investigation of the sample to sample fluctuations of thermodynamic 
quantities is nowadays feasible (whereas previously only calculation of the ensemble aver- 
age, which is less demanding computationally, was feasible). 

The particular model which is used here to study the question of the self averaging of 
fluctuations due to disorder at criticality is the bond-disordered Ashkin- Teller model on a 
square lattice. Actually this work is based on further analysis of results which were obtained 
in a previous MC study O] which aimed to determine the universality class of the model. 
The random-bond Ashkin- Teller model is particularly suitable for studying the effects of 
disorder on critical behavior. This is because the pure model possesses a line of critical 
points along which critical exponents vary continuously. In particular, the scaling exponent 
corresponding to randomness = (a/u) pur e varies continuously and is positive. Thus, 



according to the Harris criterion [12|], randomness is a relevant operator of varying strength, 
and the critical behaviour of the disordered model was indeed found to differ from that of the 
pure system. Our conclusion in the present work is that the effective exponent ratio ajv of 
a disordered model plays a central role in determining the self averaging of the fluctuations 
due to disorder at criticality. For the susceptibility, for instance, our results agree very well 
with a finite size scaling theory which we develop, according to which the relative variance 
of the susceptibility, V^-/[x] 2 , scales as l a ' v at the critical temperature. This implies lack of 
self averaging when a = (as is found for the random bond Ashkin- Teller model) and only 
weak self averaging for negative a. Our theory is successful also in describing, for models 
with weak disorder, the effect of crossover on the variance. 

Our finite size scaling theory is very similar to the physical arguments that lead to the Harris 



criterion [12], which was derived near the pure system fixed point. The difference is that we 
are assuming that similar considerations are valid near the disordered fixed point as well. 
This work is organized as follows. In section 1 we define the random bond Ashkin- Teller 



model (RBAT) and summarize its critical properties as found in a previous study fTl ]. In sec. 



PTI] we define various variances of thermodynamic quantities in disordered systems and ex- 
plain their meaning. We explain how the 'sample to sample variance' can be estimated from 
MC results. In section 3 we display our results for several bond disordered Ashkin-Teller 
models, including the four-state Potts and Ising models. We have measured the 'sample 
to sample variance' at criticality for different lattice sizes and also for different degrees of 
disorder. We discuss some qualitative features of these results, such as the apparent lack 
of self averaging and the dependence on the amount of disorder and on the specific heat 
exponent a. In section 4 we develop a phenomenological finite size scaling theory for the 
'sample to sample variance'. In section 5 we compare the predictions of the theory with the 
numerical results. We find good agreement in the case of the susceptibility and the energy, 
while the agreement in the case of the specific heat and magnetization is more questionable. 



II. THE RANDOM BOND ASHKIN-TELLER MODEL 

The model we study is the Random-Bond AT model (RBAT) on a square lattice. On 
every site of the lattice two Ising spin variables, <7j and 7$, are placed. Denoting by < ij > 
a pair of nearest neighbor sites, the Hamiltonian is given by 

H = - ^ \ K hi a i a j + K i,j T i T j + •V. ; (T,7,rr, ;/ . (5) 

<i,j> 
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The positive coupling constants ifjj and Ajj are chosen according to 

J (if 1 , A 1 ) with probability \ , , 

1 iJ ' l ' J> ~ \ (if 2 , A 2 ) with probability | ' 1 j 



The homogeneous model []13| [ (if 1 , A 1 ) = (if 2 , A 2 ) ] possesses a line of critical points, 
along which critical exponents vary continuously. This critical line interpolates between 
the Ising and four state Potts models. Even though the scaling exponent corresponding to 
randomness, (ft = (a/v), also varies continuously along this line, it takes positive values, (1 > 
> 0), so that randomness is relevant. Indeed the critical behaviour of the disordered model 



was found to be different from that of the pure one |Tl| . In JTT]] A duality transformation was 
used to locate a critical plane of the disordered model; The random model is critical when 
(if 2 , A 2 ) are the dual couplings of (if 1 , A 1 ) 0-0. This critical pi ane corresponds to the 
line of critical points of the pure model, along which critical exponents vary continuously. A 
finite size scaling study was performed for several critical models, extrapolating between the 
critical bond-disordered Ising and four state Potts models. The critical behaviour of each 
disordered model was compared with the critical behaviour of an anisotropic Ashkin- Teller 



model which was used as a reference pure model ||18|| . Whereas we found no essential change 
induced by randomness in the order parameters' critical exponents, the divergence of the 
specific heat C did change dramatically. Our results favor a logarithmic type divergence 
at T c , C ~ log/ for the entire critical manifold of the random bond Ashkin- Teller model, 
including the four state Potts model, but excluding the random bond Ising model, for which 
C ~ log log I was obtained. 

Here we give some of the details of the simulations and our main numerical results for the 
critical behavior. These are necessary for understanding and analyzing our variance results. 
All the results listed here were presented in detail in ; some essential points are reviewed 
here for completeness sake. 

Two series of critical RBAT models were studied in order to monitor two effects. The 
first series of measurements were performed at five models (or points in the couplings space), 
{(if 1 , A 1 ), (if 2 , A 2 )}j i = 0...4, which we label as Cj, i = 0...4. These were chosen so 
as to interpolate between Co, which is a random-bond Ising critical point ( A 1 = A 2 = 0), 
and C4, which is a random-bond four-state Potts critical point (A 1 = if 1 , A 2 = if 2 ). The 
points Cj interpolate in a similar manner to the way in which the critical line of the pure 
AT connects the pure Ising critical point with the pure four-state Potts critical point. The 
extent of deviation from pure behavior is determined by the difference between the two sets 
of couplings. For the series Cj, i — . . . 4 the ratio of ^ was chosen, i.e. 

if 2 « —if 1 , (7) 
10 V ; 



so that randomness will be pronounced p!9| , |20| . 

Two additional measurement points (or models) were intended to monitor the effect of 
varying the amount of randomness on the critical behaviour. The points B 2 , C 2 represent 
three RBAT models with coupling ratios ~ \ and |tt ~ |, |, ^ respectively. Thus the 
model A 2 possesses the lowest degree of randomness, while the model C 2 possesses the highest 
degree of randomness. The usual definitions for energy E, specific heat C, magnetization M 
nj, susceptibility x, polarization p = (err), and susceptibility of the polarization x^' were 
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used. Since the specific heat seems to play a dominant role in the behavior of the variance, 
we elaborate on the specific heat results, and even reproduce one graph. For the specific 
heat we found excellent agreement with the finite size scaling form 



C = ao + 6oln[l + c (Z (Q/ ^ uro -l)] 



(8) 



where (a/z/) pure is the critical exponent ratio of the corresponding Anisotropic (pure |jl8| ) 
model. Eq. (^p reproduces expected scaling forms in various limits as we now show. The 
constant Co can be expressed as 



Co 



1; 



(9) 



where l c is a crossover length, at which crossover from the pure model's power law behaviour 
to the random logarithmic behaviour occurs. Thus for I <C l c eq. reduces to the pure 
model behavior, 



C = a Q + ^(^W* , 
while for I l c and /( a / v We >la logarithmic behaviour is attained, 

C = a + b In / . 



(10) 



(11) 



Apart from crossing over to the correct pure result (^) when cq — > 0, in the Ising model 
limit, (a jv) pure — >■ 0, eq. @ becomes 



C = a + 61n(l + glial) 



(12) 



with g = co(a/V) pure . This is the finite size scaling form which was predicted analytically 
[22 1 and confirmed numerically [JT9 1 for the random bond Ising model. 
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FIG. 1. Specific heat, C, as a function of / on a log-log scale for seven critical RBAT models. 
Co is a random bond Ising model and C4 is a random bond four state Potts model. The curves are 
fits to the form ( ||) , yielding estimates for the coefficients of ( ||) which are listed in Table | . 
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In fig. D the specific heat of the critical RBAT models is plotted on a log-log scale, with 
fits to (H) using the full lattice size range 4 < / < 256. The fitting parameters a , b , and c 
together with (a/V) pure and the crossover lengths l c are listed in Table |. Note that (a/V) pure 
was not a fitting parameter, and was taken for each RBAT model from results of independent 
simulations of the corresponding anisotropic AT model. For the models Co. .4 (with large 
randomness, ^ ~ ^ ), the crossover lengths l c were found to be 1. Nonetheless these models 
differ by exhibiting some crossover with different values of (a/v) puTe (see Table |). On the 
other hand, («/V) pure of the three models A 2 , B 2 , C 2 is very similar ( .40, .37, .37 respectively) 
but they differ in their amount of randomness, -fr ~ |, \, jq respectively. Consequently, as 
one would expect, we found that their crossover lengths decrease as randomness increases: 
l c = 51 ± 7, l c = 4. ± .4, l c = 1 respectively. 

TABLE I. The fitting parameters of the critical specific heat of the random-bond Ashkin- Teller 
model, a®, bo, and Co were obtained by fitting the specific heat results of the seven critical RBAT 
models Co. a and ^2,-62 to Eq. ( ||) using lattice sizes 4 < I < 256. (^) pU re is the specific heat 
exponent of the corresponding anisotropic (pure) models. l c is the crossover length defined in (|^). 



Errors are | 


^iven in parentheses 


only when error 


is smaller than or 


of the same order as 


the number 


itself. 














a 


bo 


CO 


k 


(^)pure 


C (Ising) 


-.37(12) 


.58(1) 


5.2E4(1.5E4) 


1. 


.0001(150) 


Ci 


-4.6 


.51(2) 


1.5E6 


1. 


.171(5) 


c 2 


-4.1 


.46(127) 


5.5E4 


1. 


.375(5) 


c 3 


-3.9 


.43(4) 


5.5E4 


1. 


.549(8) 


C 4 (Potts) 


-4.1 


.42(1) 


1.0E5 


1. 


.630(8) 


B 2 


-.09(5) 


2.00(4) 


1.47(10) 


4.0(4) 


.371(5) 


A 2 


-.07(6) 


9.35(33) 


.26(2) 


51(7) 


.40(1) 



We found that the magnetization M, susceptibility x an d the susceptibility of the po- 
larization, x are wei l described at criticality by the following scaling laws: 

M = A M r» , (13a) 
X = A x li , (13b) 

X ip) = A x(p) l^- . (13c) 

The estimates for the exponents ^, ^ and which were obtained using lattice sizes / > 24, 
are listed in Table ||. Even though one observes slight variation of ^ and ^ from model to 
model, the results are consistent also with fixed, non- varying exponents f = | and ^ = |, 
modified by a logarithmic correction. So ^ shows very little variation or does not vary at 
all. This is nearly the same behavior as was found for the corresponding anisotropic models 
where 2 is predicted analytically [£23, [23] to be constant for all models ^ = ?. The exponent 



(p) 

ratio ' J -^- connected with the susceptibility of the polarization which varies continuously for 
the pure Ashkin- Teller model seems to do so also for the random models (see Table !]]). 
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TABLE II. Critical exponents ratios for seven critical RBAT models Co. .4 and A 2 ^Bi- These 
exponent ratios for the magnetization M, the susceptibility \ an d the susceptibility of the polar- 
ization x were obtained by fitting results for lattice sizes I > 24 to equation ( p^). 





1 

V 


£ 

V 


7 (p) 

V 


C (Ising) 


1.751(5) 


.125(3) 


1.549(9) 


Ci 


1.751(6) 


.124(3) 


1.575(8) 


c 2 


1.743(5) 


.129(3) 


1.597(9) 


c 3 


1.736(3) 


.133(2) 


1.638(5) 


C 4 (Potts) 


1.714(5) 


.145(3) 


1.714(5) 


B 2 


1.738(4) 


.132(3) 


1.586(6) 


A 2 


1.739(5) 


.132(3) 


1.590(8) 



III. VARIANCES; DEFINITIONS AND ESTIMATORS 

In this section we define two types of variances of thermodynamic quantities in disordered 
systems and explain their relation to error estimates. We explain how the 'sample to sample 
variance' can be estimated from MC data. 

First consider some sample x which is simulated at some temperature T. Because of the 
thermal fluctuations and the finite simulation time, we obtain for this sample an estimate 
P x of the exact P x , with an error 

2 

{8P~ X ) 2 = -^-, T MC large. (14) 

J-MCHx 

Tmc is the length of the MC runs and r is the autocorrelation time of the MC dynamics. 
a\ x is the variance of P within the sample x due to thermal fluctuations. In practice, in 
order to avoid the estimation of r which requires a long simulation time, we estimate (SP X ) 2 
by binning the MC sequence into ~ 10 subsequences and treating each subsequence as 
independent ( the Jack-knife procedure). 

The estimate for the error in the estimation of [P x ] , the average of P over all samples, is 
given by 

1 n 

(m]) 2 = 7 TT-E(^-[P,]) 2 , n large, (15) 



{n — l)n 



x=l 



where n is the number of random bond samples. In contrast to fll4"D this total error has two 
contributions, namely the sample to sample fluctuations of the exact P x around [P x ] and 
the thermal fluctuations of P x around P x within each sample,that is 

(S[P x ]) 2 = - + [ ° T > X ], n,T MC large. (16) 

n uTmc/t 



Thus by estimating (5P X ) for all x and (5[P X ]) with (fL5D we obtain V through (16); it is 



an unbiased estimate of the variance of the exact P x due to sample to sample fluctuations 



S 



(see ref. p5| for a basic statistical explanation). In order to minimize the error of [P x ] for 
a given amount of computer time, one needs to adjust T MC so that the two terms in flTEp 
are equal. However, if one is interested in obtaining a reasonable estimate of V, T MC needs 
to be chosen larger, so as to obtain accurate estimates of the P x 's and minimize the second 
term on the l.h.s. of (|T6|). 

As explained in the introduction, the dependence of thermal variance <y\ on the lattice 
size I has been examined (for homogeneous models) in ref. 0,0. Thus from here on the 
term variance will refer to the variance due to disorder. Here it is our aim to examine the 
dependence of the variance V on I at criticality, one reason being that for MC simulations of 
disordered systems, it has the bigger influence on their accuracy. This is in addition to the 
theoretical motivation given in the introduction. In the next section we display our results 
for the variance V of the random-bond Ashkin- Teller model. 



IV. VARIANCE RESULTS OF THE RANDOM BOND ASHKIN-TELLER 

MODELS 

A. Far from Criticality 

Far from criticality the correlation length is finite and one would expect the system to 
behave similarly to a collection of independent smaller systems. Thus one would expect the 
Brout argument to hold with the variance scaling as l~ d . Nonetheless this is not obvious: 
Note that the thermal fluctuations of the specific heat C and the susceptibility \ are non 
self averaging even away from criticality RID]. Thus the RBAT model C2 (the choice 



of model was arbitrary ) was simulated at the reduced temperature t = 1. . In Fig. |2| 
we show the relative variances Vp/[P} 2 , where Vp is defined in (|16|) and P = E, M,C,x, 
as a function of log 10 Z. The linear curves are fits to the form Vp/P 2 = Al~~ p . We find 
p = 2.06(7), 2.13(7), 2.04(6), 2.12(7) for X ,C,E,M respectively. Thus the Brout argument 
is confirmed and far from criticality strong self averaging holds. 
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FIG. 2. The relative variances V P /[P] 2 , P = E, M, C, x as a function of log / for the RBAT 
model C2 at the reduced temperature t = 1. The linear curves are fits to the form Vp/P 2 = Al~ p , 
yielding p = 2.06(7), 2.13(7), 2.04(6), 2.12(7) for X , C, E, M respectively. 



B. The Variance at Criticality 

1. Distributions 

In order to visualize how large the sample to sample fluctuations are, at the critical 
temperature, several histograms of the frequency of occurrence of samples according to their 
susceptibility Xx or according to their specific heat C x , are shown in figures ^|-[F[ The abscissa 
is scaled by the average susceptibility [x] (or specific heat [C]) of all samples. The histogram 
of the susceptibility for lattice size I = 192 is shown in Fig. |3] for the Ising model Co and 
in Fig. [| for the four-state Potts model C4. The frequency scale of both figures is scaled so 
that the area of both histograms is the same. Even though the lattice size is rather large, 
the distributions are very wide; a measurement of a value of \ at 40% above the mean [%] 
has a non-negligible probability for the four state Potts model. 
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FIG. 3. Histogram of the frequency of occurrence of samples according to their susceptibility 
scaled by the average susceptibility, for the Ising model Co and lattice size I = 192; with 240 
samples. 



O . 5 1 1.5 2 

FIG. 4. Histogram of the frequency of occurrence of samples according to their susceptibility 
scaled by the average susceptibility, for the four-state Potts model C4 and lattice size I = 192; with 
370 samples. 
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FIG. 5. Histogram of the frequency of occurrence of samples according to their susceptibility 
scaled by the average susceptibility, for the four-state Potts model C4 and lattice size I = 24; with 
920 samples. 

There is a marked difference between the width of the distribution of the Ising model 
{JVx/\x\ ~ 0.2 ) and the much wider distribution of the four-state Potts model {Jv^/\x\ ~ 
0.32). The histogram of the susceptibility for the four-state Potts model with lattice size 
I = 24 is shown in Fig. [5[ Note that the width of the distribution here is slightly narrower 
(« 0.29) than that of Fig. |. This very small difference (and even slight increase) of the 
width as I increases hints at a lack of self averaging of the susceptibility of the four-state 
Potts model. An additional striking difference between the susceptibility distributions of 
the four state Potts model, figures |] and |5|, and the Ising model, Fig. |3], is that the former 
are strongly asymmetric (this asymmetry was measured by measuring the third moment of 
the distribution). A possible explanation for this asymmetry, which exists to some degree 

[V] and in 



in all the models, is given in Sec. 
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The average errors in the estimation 
of the susceptibility of a single sample x, divided by the average susceptibility, [^Xs]/[x], 
are 0.01, 0.017, 0.012 in figures |3], [| and || respectively. Since these errors are negligible as 
compared to the widths, the histograms are highly reliable. 

The histograms of the specific heat for lattice size I = 48 are shown in Fig. ^| for the 
Ising model Co and in Fig. [7| for the four-state Potts model C4. Note that the distributions 
of the specific heat are much narrower than those of the susceptibility. The width of the 
distribution for the four-state Potts model (~ 0.126) is about twice wider than the width of 
the distribution for the Ising model (~ .062). The asymmetry of the distribution for the four- 
state Potts model is almost unnoticeable and is of the opposite sign than the asymmetry 
of the susceptibility. The average error in the estimation of the specific heat of a single 
sample x divided by the average specific heat [5C^]/[C] is 0.038 for the Ising model and 
0.048 for the four-state Potts model, so that these histograms are much less accurate than 
those of the susceptibility. For the larger lattices the ratio between the width and the error 
becomes smaller, mostly because the width becomes smaller, and histograms become even 
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less accurate. Thus in order to obtain more accurate histograms and also better estimates 
of the variance (which is the square of the width of the histograms), longer simulation times 
would be needed, in order to obtain more accurate estimates of the C x 's. This may be done 
in a future study. 
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FIG. 6. Histogram of the frequency of occurrence of samples according to their specific heat 
scaled by the average specific heat, for the Ising model Co and lattice size / = 48 with 600 samples. 
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FIG. 7. Histogram of the frequency of occurrence of samples according to their specific heat 
scaled by the average specific heat, for the four-state Potts model C4 and lattice size / = 48 with 
630 samples. 
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2. The variance 



In Fig. [S] we show the variance of \i °f the seven critical RBAT models. For the sake 
of clarity (so that the data do not fall on top of each-other) V x of the model was multiplied 
by 2 i+1 . The lines are fits according to a theory which we develop in the next section. Here 
we just note that V x is measured with high precision, so that it may be faithfully tested 
against theory. 
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FIG. 8. The variance of x> V x as a function of log / for all critical models, Co. .4, and A 2 ,B 2 
of the RBAT model. For the sake of clarity V x of the model Cj was multiplied by 2 l+1 . The solid 
lines are fits to the form ( [|l|) , yielding estimates for fitting parameters which are listed in table 

The relative variance V x /[x} 2 is plotted in Fig. |9|. Since it is the ratio of two fluctuating 
quantities, the errors are quite large. Nonetheless the main trends can be seen. First note 
that apparently for all models (except for the weakly random model A 2 ) V x /[x] 2 — > const, 
so that the susceptibility is non-self averaging. It is also possible that V x /[x} 2 is slightly 
increasing with / for some models ( e.g. the four state Potts model C4 ) or slightly decreasing 
for the Ising model Co- Upon comparison of the models Cj i = . . . 4 we make the following 
observations. The higher is the specific heat of a model, the larger is its relative variance 
(see Fig. [I]). The higher is the exponent (-) pure of the pure model (see Table |), the larger 
is the initial slope of the relative variance of the corresponding random model. Thus the 
relative variance of the Ising model Cq is the smallest and the increase with I, for small I, 
is the smallest. The relative variance of the four-state Potts model C4 is the largest and 
the increase with I, for small I, is the largest. The relative variances of the RBAT models 
Ci,2,3 fall in between. The relative variance of the weakly random model A 2 shows a steady 
increase with Z, in contrast with the highly random model C 2 , in which a shorter increase 
is followed by a plateau. This is reminiscent of the specific heat of the A 2 model which 
exhibits very slow crossover from the power-law behavior ( |T0D to the asymptotic logarithmic 
behavior ( p|) with a crossover length of l c ~ 50. Thus for small lattice sizes the model 
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A2 exhibits effective exponents (of the specific heat) of the pure model, and also exhibits 
a small variance due to its small degree of randomness. But as the lattice size increases, 
this effect diminishes, the effective exponents approach the random value, and the variance 
approaches that of the highly random models. 
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FIG. 9. The scaled variance of the susceptibility, V x /[x] 2 as a function of logZ for all critical 
models, C0..4, and ^2,-62 of the RBAT model. 

A very similar picture is obtained for the relative variance of the magnetization Vm/[M] 2 , 
as seen in figure 10]. The qualitative picture of the magnetization results, Fig. |1(| is very 
similar to that of the susceptibility results, Fig. ||, showing the same trends as outlined 
above. Yet we emphasize that even though the magnetization is an intensive quantity, it 
does not seem that Vm/[M} 2 — >• as I increases so that the magnetization is not self 
averaging at criticality. 
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FIG. 10. The scaled variance of M, Vm/[M] 2 as a function of log / for all critical models, Co. .4, 
and A 2 , B 2 of the RBAT model. 



In Fig. |TT| the variance of the energy Ve |27| ^ s plotted on a log-log scale. For the 
sake of clarity (so that the data do not fall on top of each-other) Ve of the model Cf was 
multiplied by 2 l+l . We fit the data to the form Ve ~ l~ e for lattice sizes I > 16 (but the 
fitting curves shown in Fig. |ll| are not made with this form but with a more complicated 
one which is due to a theory which we develop in the next section). The highest value of 
9,9 = 1.855(13), was obtained for the Ising model C . For the four-state Potts model, C 4 , 
we obtained 9 = 1.72(2), while for the models Ci^,3 the values of 9 fell between these two 
values. Thus in contrast with the susceptibility and the magnetization, the variance of the 
energy Ve is weakly self-averaging. But similar to the susceptibility, models with a higher 
specific heat or with a higher effective - have a smaller 9, and thus their Ve decreases more 
slowly with I. For the weakly random model A 2 , 9 = 1.30(3) so that again the slope of the 
variance is correlated with the high slope of the specific heat of this model. 
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FIG. 11. The variance of the energy, Ve as a function of log? for all critical models, Co.. 4, and 
A2, B2 of the RBAT model. For the sake of clarity Ve of the model was multiplied by 2* +1 . The 
solid lines are fits to the form ( ^) , yielding estimates for the fitting parameters a v , b v which are 



listed in table IV . 



The results of the relative variance Vc/C 2 , plotted in Fig. 12, seem to indicate that the 
specific heat is weakly self averaging. Nonetheless the effective slopes increase with I (or the 
absolute values of the slopes decrease with I, this trend being strongest for the four-state 
Potts model C4 ) so that it is possible that self averaging does not hold for very large I. It also 
seems possible that the Ising model is self-averaging while the other models are not. Clearly 
more accurate data and data from larger systems would be useful. As in other variances, we 
observe qualitatively that the relative variance of the moderately random models, A 2 and 
B 2 , approaches that of the highly random ones as I increases and even exceeds it. 
The findings of this section are partly summarized in table [[II], where the self averaging 
properties of the highly random critical models are displayed. 
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FIG. 12. The scaled variance of C, Vc/[C] 2 as a function of log£ for all critical models, Co. .4, 
and A 2 , B 2 of the RBAT model. 



TABLE III. Summary of the self averaging properties of the critical random-bond four-state 
Potts C4, Ashkin- Teller C1..3, and Ising Co models. The letter 'n' stands for non self averaging, 'w' 
for weakly self averaging, and '?' for inconclusive results. This summary is according to a subjective 
examination of the numerical results as displayed in figures P , [To| , pj] , |l2 13 . According to our 



theory only the energy E is weakly self averaging while all other quantities are non self averaging 
in all of these models. If the theory is correct then other behavior implied by the numerical data 
are merely transients. 



model 


X 


X (P) 


M 


C 


E 


C (Ising) 


? 


? 


? 


w 


w 


C1..3 (Ashkin-Teller) 


n 


n 


n 


w 


w 


C 4 (4 state Potts) 


11 


n 


n 


w 


w 



In the next section we develop a phenomenological finite size scaling theory for the 
variance. This theory explains the apparent connection between the variance and the specific 
heat behavior of the random models. In the last section we explain how this theory was 
applied to the results we have displayed here and discuss the comparison between our scaling 
theory and the numerical results. 



V. FINITE SIZE SCALING OF SAMPLE TO SAMPLE FLUCTUATIONS AT 

CRITIC ALITY 

As our numerical results show, we have obtained quite accurate estimates of the variance 
V of the thermodynamic functions at the critical temperature for different lattice sizes I. 
In order to understand these results, a phenomenological theory of finite size scaling of 
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disordered systems, which will take into account sample to sample fluctuations, needs to be 
developed. 

The main result of our theory will be the scaling of the variance V with I at criticality. 
To be precise, we will calculate the variance V of P (e.g. P = C, M, E or x, where all these 
quantities are normalized per volume; i.e. they are densities) 

V(T,l) = [(P x (T,l)-[P x (T,l)]) 2 \- (17) 

P x (T,l) is the exact value of P (that is, after the thermal fluctuations have been averaged 
over) of a specific sample x (with some specific realization of randomly distributed bonds) of 
linear size I at temperature T. Again the square brackets denote averaging over the different 
samples x. 

Our conclusion will be that when the specific heat exponent a is negative the leading 
behavior of V at T c is 

V(T e ,l)~K* r l«+!, implying ~ • (18) 

Where K r is a measure of the amount of randomness or disorder and p is the critical exponent 
of the quantity P, e.g. if P — % then p = ' 1 . Eq. (|i8|) implies that disordered systems at 



criticality are only weakly self averaging when ^ < 0. For ^ = 0+ (log), as was found [TT 
for the random bond Ashkin- Teller model, our derivation is strictly not valid for I > 1. 
Nonetheless for the range of lattice sizes considered, we found good agreement between 
the numerical results for the variance of x, Xp an d E and theoretical fits according to (|T 



together with next to leading terms (see figures [8] , and |ll] and discussion in the next 
section). If no dramatic change occurs at larger sizes, then the sample to sample fluctuations 
of the random bond Ashkin- Teller model are non self averaging. 

The result ( |TE| ) indicates that the sample to sample fluctuations at the critical temper- 
ature T c depend strongly on the specific heat exponent -. This strong dependence can 
be made plausible based on heuristic arguments. These heuristic arguments will serve to 
define some basic ingredients of our approach and will be followed by a more quantitative 
treatment. 

We start by characterizing every specific sample x of size I by a pseudo-critical temper- 
ature T c (x,l). This pseudo-critical temperature can, for instance, be the temperature at 
which a maximum in the specific heat of the sample occurs. We denote the average pseudo- 
critical temperature as T C {1) = [T c (x, I)]. We assume that, as is the case in homogeneous 
systems, 

T c {l)-T c = ar yt , (19) 

where a is a constant, y t = 1/V and T c = lim/^oo T c (/). T c is the average critical temperature 
of the ensemble of infinite samples. Eq. ( |T§D is supported by a numerical study of the 
three dimensional dilute Ising model. 

Next we assume that T c (x, I) fluctuates around T c (l) with width 

5T c (l) ~ l~ d/2 . (20) 



This assumption is probably true [j|,[nj for small disorder and small I, or close to the pure 
system fixed point. We assume it without proof, for large disorder as well, though for large 
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disorder (or close to the random fixed point) the possibility that ST c (l) ~ I Vt has been 
raised ||. 

Define reduced temperatures 

t c (x,l) = Tc{xJ ^' Tc , (21a) 

J- c 



T c (l) - T c 

and the reduced width 



t c (l) = CV T C , (21b) 



St c (l) = 6 -W- ■ (21c) 



T, 



c 



We make (|20|) more specific by assuming for t c (x, I) a Gaussian probability distribution 

q(t c (x,l)) 

jd/2 

q(t c (x,l)) = __exp{-[t c (x,0 -t c (/)] 2 /72^ r 2 } • (22) 

V ZlXK r 



The width of the distribution is controlled by the lattice size / and by K T which is a measure 
of the amount of randomness or disorder. 

The scaling relations |20D already make the result fll8|) plausible. The main idea 
is that the sample to sample fluctuations at T c are governed by the relative magnitude of 
two temperature differences. The first is the difference between the average pseudo-critical 
temperature T C {1) and the critical temperature of the infinite system T c . The second is the 
difference between T c (l) and T c (x, I), the pseudo-critical temperature of the sample x, which 
is governed by St c (l). If 5t c {l) ^> \t c (l)\ then fluctuations in t c (x,l) are so large that for 
some samples one will find T c > T c (x, I) while for other samples T c < T c (x, I). In this case, 
even though we are simulating all samples at T c , some samples are in their high temperature 
phase while others are in their low temperature phase. This will obviously increase the 
sample to sample fluctuations in any observable. If, on the other hand, St c (l) <C |£ C (0I> then 
T c — T c (x, I) will always have the same sign and fluctuations will be smaller. The condition 
8t c (l) <C \tc{l)\ w iH be fulfilled for large / if y t — f < or, using the hyper-scaling relation 
- = 2y t — d, if - < 0. For disordered systems, the bound y t < d/2 has been proven by Chayes 



et al. ||30 |, so that asymptotically one always finds - < 0. However, for small I and small 
disorder, the system may be governed by a positive (^) pur e- In this case sample to sample 
fluctuations can increase with lattice size, as is indeed seen in our numerical results for the 
weakly disordered model A 2 . Thus on the basis of these considerations one can conclude 
that the sign and magnitude of the specific heat exponent a of the disordered model have a 
strong influence on the sample to sample fluctuations |2(J , and will determine whether they 
are self averaging. The discussion above is analogous to the physical arguments leading to 



the Harris criterion |12], but in a finite size scaling formulation. The difference is that the 
Harris criterion was derived near the pure system fixed point, while we are assuming that 
similar conditions apply also next to the disordered critical fixed point. 
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In order to put these general considerations on more quantitative grounds, we proceed 
to derive the finite-size scaling expression (|18|) for the variance of various thermodynamic 
quantities. Start by introducing the reduced temperature of each sample x; 

i. = ■ (23) 

J- c 

We assume ( for samples with T close to T c (x, I)) a finite size scaling form for the singular 
part of P x , 

Pr s {T,l) = l»Q x {i x V*) . (24) 

The form of the function Q X (Z) (or its coefficients) are assumed to be sample dependent 
but the critical exponents p, y% are assumed to be universal or sample independent. 



eq. (p4|) embodies the usual |3T[ finite size scaling assumption that in the vicinity of the 
critical temperature the behaviour of a large finite system is governed by the scaled variable 
£x/l- We use this assumption, even though in the present context it implies that a single 
correlation length £ x is sufficient to describe the state of a disordered sample, which is not 
obvious at all. This "thermal" I -dependence is compounded by the fact that if we increase I, 
we must generate additional random bonds, and hence increasing I necessitates, effectively, 
changing x (that represents a particular realization of the random bond variables). Since 
x affects P^ in 9 through the non- universal coefficients of Q x , a non-thermal dependence of 
psmg on £ - g i nc i ucec [ m ain task of our analysis is to separate the thermal I dependence 
from the non-thermal component. 

At this stage it is possible to draw some more conclusions based on (p4]), without making 
strong assumptions about the coefficients of Q x . We leave such derivations for the Appendix. 
Here we proceed in a more straightforward manner by using a simplifying ansatz. Our ansatz 
states that the coefficients of Q x depend only on At c (x, I), the deviation of the pseudo-critical 
temperature of the sample from the average pseudo-critical temperature, defined as 

At c (x,l) =t c (x,l) -t c (l) . (25) 

It is convenient to proceed by rewriting i x as i x = t — At c (x, I) —t c (l) with t = T ~ Tc . Using 
the scaling of t c (l) [see ( |19|) and (|21b|) 1 , we substitute Q x by a different scaling function Q x 



and rewrite Eq. (^) as 

psin 9{Ti t) = lP Q x{{f _ Atc(X) /} _ tc (/)) P } = lPQ x{ ( t _ Afc ( Xj . (26) 

For completeness of the treatment which will later prove to be necessary we do not neglect 
the analytic dependence of P X (T, I) on (t — At c (x, I)) JHJ, and write 

P X (T, l) = A x + B x (t - At c (x, I)) + C x (t - At c (x, I)) 2 + ... + l p Q x {(t - At c (x, /))P} . (27) 

The coefficients A X ,B X ,C X are assumed to be sample dependent in the same way that the 
coefficients of Q x are; namely they depend only on At c (x,l) ]33|]. Next, assume the depen- 



dence of the coefficients on At c (x,l) is analytic. Since according to ( |22|) and (p5|), At c (x,l) 
is distributed around zero with width that scales as /~ d / 2 , we can expand 
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A X = A + A t At c {x, I) + A 2 (At c (x, I)) 2 + ... , (28) 

where Ao,A±, A 2 are sample independent. The same type of expansion is assumed for B x , C x 
etc. 

We are interested in knowing what happens at T = T c , the average critical temperature 
of the ensemble of infinite samples. Thus we set T = T c which implies t = 0. For the analytic 
part of fl27D we get 



p analytic (Tc; q = ^ _ /} + Cx At c (x, I) 2 + . . . = 

A + (A, - B )At e (x,l) + (A 2 -B, + C )(At c (x,l)) 2 + ... 

= a + bAt c (x,l) + c(At c (x,l)) 2 , (29) 

where the second equality is reached by use of (|28| ) and the same expansions for other 
coefficients. The last equality is a redefinition of constants. In a similar way we expand Q x 

Q X {Z) = D X + E X Z + F X Z 2 + ... , (30) 

where D X ,E X ,F X are again expanded as in ( pS|) . Again setting t = 0, we obtain for the 
singular part of ( p7f) 

psingular^ q = + ^ _ ^p^^ /) + (£> a _ E JVt + F l 2yt )(At c (x, I)) 2 } + ... . 

(31) 

We stress that since we set Z = —At c (x, l)l yt , Z is fluctuating around zero with width that 
scales as ~ . Thus the expansion ( p0|) is justified asymptotically only for a < 0. Putting 
together (p9|) and (|3lf) we have 



P X {T C , = (a + £W P ) + (6 + ^ p - ^ / p+yt )^c(x, + 

(c + D 2 Z' - E x \p^ + F /^)(At c (x, Z)) 2 + . . . = d + eAt c (x, I) + /(At c (x, I)) 2 • (32) 

Notice that here the only dependence on the specific sample x is through explicit depen- 
dence on At c (x, I), the deviation of its reduced pseudo-critical temperature from the average 
pseudo-critical temperature. Taking the quenched sample average [ ] with the probability 
distribution (221), using [At c (x, I)] = 0, we get 



[P x (T c ,l)]=d + f[(At c (x,l)) 2 ], (33) 
and using [(At c (x, I)) 3 } = we further obtain 

[(P,(T C ,/)) 2 ] = d 2 + e 2 [(At c (xJ)) 2 } + f 2 [(At c (x,l)) i } + 2df[(At c (x,l)) 2 } . (34) 
The variance is then given by 

V(T C , I) = e 2 [(At c (x, I)) 2 ] + f 2 {[(At c (x, /)) 4 ] - [(At c (x, I)) 2 ] 2 } = e 2 [(At c (x, I)) 2 ] + 2f 2 [(At c (x, I)) 2 ] 2 

(35) 

where the last equality is a property of the Gaussian distribution. Lastly we use 
[(At c (x, I)) 2 } = K 2 l~ d and obtain to the leading orders in I 
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V{T C , I) = (b 2 + D 2 l 2p + E$l 2p+2yt + 2bD 1 F - 2bE F +yt - 2D 1 E l 2p+yt )K 2 l- d + 

( F 2fp+*vt)2K?r 2d + . . . . (36) 

Since y t > 0, and usually p + y t > 0, the leading term in (^) is 

V(T C , I) ~ E 2 ^ 2 / 2 ^ 2 ^ = E 2 K 2 l 2p+ ^ . (37) 

The last term in fl36|) is proportional to K^l 2p+2 ^ , and may be neglected with respect to 
(|57|) only for ~ < 0, or if K 2 <C 1 and I is not too large. ([37]) is our main result for the 
variance, where all exponents p,a,i> are exponents that characterize the disordered system. 
It means that disordered systems at criticality are only weakly self averaging when - < 0. 
Though our derivation is not valid when ^ = 0+ (C ~ log/) it seems that in this case there 
is no self averaging of the sample to sample fluctuations. This point is further discussed in 
the next section. 

We note that for a < and in the large I limit considered at the end of the Appendix 
[equation ( |A4T ) ], where 

Q x (-At c (x, Z)P) -> Q(-At c (x, /)P) , (38) 

the coefficients of Q x are independent of x so that D x = D , E x = E etc. . Neglecting the 
analytic part of P, this limit corresponds to our derivation with only D , E , F ^ and all 
other coefficients ( Di,D2,E\ etc. ) equal to zero. Thus in this limit the main result (|37D 
is unchanged, though less assumptions are needed. 

From ([33]) corrections to the scaling of P are obtained, 

[P X (T C , l)] = a + D l p + (c + D 2 l p - E x l p+yt + F l p+2yt )K 2 r d ■ (39) 

So that the leading behaviour of P is 

[P X {T C , l)]=a + D l p + F K 2 r l p+2y t- d = a + D F + F K 2 F + ^ . (40) 

Thus for negative a, the third term in ( |40"D is a correction to scaling due to sample to sample 
fluctuations. It follows that for ^ < 0, (|0]) and (|37|) are consistent with V/[P] 2 ~ 1*. A 
special case is when a pure < and randomness is an irrelevant operator (at the pure system 
fixed point) with a scaling exponent (^) pU re- In this case the disordered system has the 
same exponents as the pure one, ^ = (^) pU re- Therefore the correction to scaling we have 
obtained due to sample to sample fluctuations has the same exponent as the correction term 
connected with the irrelevant operator corresponding to randomness. 



VI. COMPARISON OF THEORY WITH VARIANCE RESULTS OF THE RBAT 

MODELS 

The derivation presented in the previous section as can be readily seen from equations (|27| , 



28 , |30|) , involved an expansion in the two parameters J[(At c (x,l)) 2 ] and J[(At c (x, l)) 2 ]l yt 



These scale as K r l~i and K r F f ~2 = K r l^ . Thus the derivation is valid for small K r , 
meaning small disorder and small -. For negative ^ the validity of the expansion improves 
as I increases, while a positive - is not possible |3(| . 
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In the case of the random bond Ashkin Teller model, we have asymptotically C ~ b log I 
so that - = + . It seems that in this case the expansion is not justified. Practically though, 
for the accessible range of lattice sizes, things depend on the constant of proportionality b. 
If b is small, then for a finite but large interval of lattice sizes / the expansion is justified. 
Indeed, in the case of the highly random RBAT models Co.. .4, b falls in the range 0.138(4) < 
b < 0.280(6). Upon inspection of Fig. one may also see that the value of the specific heat 
of these models shows very little variation for lattice sizes I > 16. Thus the parameter K r l~ , 
which should scale with I as the specific heat does, increases very slowly with I. This implies 
that for the accessible range of lattice sizes I our expansion is valid. The specific heat of the 
weakly random model A 2 effectively diverges with a positive effective ^ but because of its 
weak degree of randomness there is good reason to believe that the expansion will be valid 
due to a small value of K r expected for models with small randomness. The B 2 model with 
moderate disorder is expected to fall between the A 2 model and the Co.. .4 models. Thus 
there is reason to hope that our theory is applicable to the variance results in the accessible 
range 4 < I < 256. Indeed the agreement we now display between numerical data and theory 
is good. 

For observables with p > the two leading terms in ( |3~6"D are the third and sixth terms. 
We use hyper-scaling to write l yt 
specific heat (T 



1 2^ + 2 and substitute in (|T8|) l» by the behaviour of the 



V(T c ,l) 



a 



Thus we propose for the RBAT models the leading behaviour |34 
u Z 2 "ln[l + c (Z {c ^ )puro - 1)] + c„Z 2 ^{ln[l + c (/ (a/iy)purc - 1)]}* , 



(41) 



with a v = EqK% and c v = —2D\EqK^. (note that for every thermodynamic quantity there 
are different coefficients Eq,D\ etc. ) The expression ln[l + co(Z <Q//i/ - )puro — 1)] describes the 
singular behavior of the specific heat including the crossover from the pure model behavior , 
characterized by the pure model exponent (a/z/) pure (see Eq. (||) )• p is the critical exponent 



of the quantity whose variance is measured (e.g. 



7 



for x 



(p) 



In fig. 



we show the variance of x, V x 



P — „ f° r X an d p 
of the seven critical RBAT models St- 
and - were taken from Ta- 



ted by the function fl4"T|), where the parameters c , (a/z/) pure 
bles I and [TT]. For the sake of clarity (so that the data do not fall on top of each- 
other) V x of the model C was multiplied by 2 l+1 . The fitting parameters 
given in Table |Iy|. The agreement with our scaling prediction is quite encouraging. 



TABLE IV. Fitting parameters for the variances of x, X \ E f° r the critical models Co. .4, and 
A2, B 2 according to eq-s. ( |4l| ) and ( f42| ), using lattice sizes I > 8. 





v x 

a v 


c v 


V x (v) 

a v 


C v 


V E 

a v 


K 


Co(Ising) 


0.0145(7) 


0.11(2) 


0.033(1) 


.07(4) 


0.29(7) 


0.8(2) 


Ci 


0.0039(1) 


0.026(10) 


0.0134(3) 


-0.03(2) 


0.128(14) 


-1.45(19) 


c 2 


0.0059(1) 


0.014(10) 


0.0172(2) 


-0.09(2) 


0.135(13) 


-1.19(16) 


c 3 


0.0069(2) 


-0.01(1) 


0.0147(3) 


-0.12(2) 


0.133(15) 


-1.27(18) 


C 4 (Potts) 


0.0082(2) 


-0.04(1) 


0.082(2) 


-0.04(1) 


0.121(13) 


-1.25(17) 


B 2 


0.033(1) 


0.13(2) 


0.092(2) 


0.19(4) 


0.99(15) 


1.2(2) 


A 2 


0.056(1) 


0.028(12) 


0.143(3) 


0.05(3) 


5.49(44) 


1.28(23) 
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The same analysis has been carried out for the variance of x^\ where was taken from 
Table [TI], and the results are plotted in fig. [13]. Again the fitting parameters a v , c v are given 
in Table |IV| and the agreement between the numerical results and our scaling prediction is 
encouraging. We stress that the only fitting parameters of the fits in figures [8] and ^ are 
a v c„; the other parameters of eq. (|4l|), p, Co, (ct/^) p were obtained previously [Tl[] from the 
specific heat results and from the results for \ and x ■ 
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FIG. 13. The variance of x^ p \ K^Cp) as a function of logl for all critical models, Co.. 4, and 
A 2 ,B 2 of the RBAT model. For the sake of clarity V x ( v ) of the model Cj was multiplied by 2 l+ . 
The solid lines are fits to the form ( ^T|) , yielding estimates for fitting parameters which are listed 



in table IV 



Since the first term in ( ffl]) is the dominant one (by a factor of l yt , where yt > 1 ), we test 
(p|) again in another manner, in fig. [14] we plot the scaled V x : V x l~ 2p / ln[l+co(/^ a ^- )puro — 1)]. 
Indeed it seems that the data points approach a constant value, confirming the leading term 
IT) which originates in the leading behavior of the variance fl3~T|) . 
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FIG. 14. The scaled variance of x> V x l 2 ^ / m [l + co(^ Q//l/ ^ p — 1)] as a function of log / for all 
critical models, Co. .4, and A2,B2 of the RBAT model. 

For the energy p = (a — 1)/V < so that the two leading terms in (|36D are the third 
and the fifth ones. Again by using hyper-scaling and substituting by the behaviour of 
the specific heat (M), we arrive at the scaling form for the variance of the energy 



V E (T C , I) = K(ln[l + c (Z (Q/ ^ c - l)]) 2 + b v ln[l + co(l {a/u)vm - 1)]}Z" 



(42) 



with b v = —2bE^K 2 r In fig. |ll] we show the variance of the energy, Ve of the seven critical 
RBAT models fitted by the function ([42]). For the sake of clarity Ve of the model Ci was 



multiplied by 2 i+1 . The agreement between theory and the numerical data is good and the 
fitting parameters a v , b v are given in table |V|. 

For the magnetization and the specific heat p — - and p = — respectively. In these 
cases \p\ is small and the fifth and the sixth terms in Eq. (|36"D are of similar order in I . 
Thus one may not neglect one term with respect to the other as was done for the energy 
and the susceptibility. Thus we fit the variance of the specific heat to the form 



V C (T C , I) = a v {\n[l + c Q (l (a/ ^ - l)]} 3 + W~5{ln[l + c (l^— - 1)]} L5 + 



c v l 2{ln[l 



c (l {a/v)p ™ 



1)]} 



2.5 



(43) 



In fig. ^ the variance of the specific heat, Vc of the seven critical RBAT mod- 
els is fitted by the function ([|3|), with the fitting coefficients given in Table [V]. 
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TABLE V. Fitting parameters for the variances of the specific heat Vc and the magnetization 
Vm for the critical models Co.. 4, and A 2 ,B 2 according to eq-s. (|43| ) and ), using lattice sizes 
I > 4. 







V C 






V M 






a v 


by 


c v 


a v 


by 


Cy 


C (Ising) 


0.000016(12) 


0.028(10) 


0.008(4) 


0.0037(1) 


0.25(2) 


-0.28(2) 




0.0000011(1) 


-0.045(4) 


0.0038(3) 


0.00108(3) 


0.102(9) 


-0.12(1) 




0.0000028(5) 


-0.078(8) 


0.0081(7) 


0.00169(7) 


0.10(2) 


-0.12(2) 


c 3 


0.0000079(6) 


-0.069(5) 


0.0069(5) 


0.00213(9) 


0.054(20) 


-0.068(24) 


C 4 (Potts) 


0.0000099(4) 


-0.062(5) 


0.0057(4) 


0.0027(1) 


-0.004(23) 


-0.004(28) 


B 2 


0.0079(6) 


-0.44(9) 


1.4(1) 


0.0089(3) 


0.27(3) 


-0.30(4) 


A 2 


0.67(7) 


-4.9(18) 


36.5(54) 


0.0167(8) 


0.07(7) 


-0.07(10) 



The data for large lattice sizes is rather noisy and three parameter fits are not so reliable 



with only eleven data points, so that both the results and the fitting curves in Fig. [L5 
should be taken with a grain of salt. The obtained fitting coefficients are consistent with the 
coefficient E being much smaller than b, D\. A small value for E is quite plausible if the 
specific heat as a function of the temperature is close to being symmetric around the critical 
point [see (|30l)]. This symmetry is supported by the symmetric form of the histograms of the 
specific heat as shown in figures |^ and [F| For the Ising model Co the errors of the coefficients 
a v , h v and c v are of the same order of magnitude as the coefficients themselves. However for 
the other models the errors are reasonable and though a v is small, we have EqK^ = a v > 0, 
meaning that, asymptotically, the first term in fl43D will dominate. This implies that the 
specific heat of the RBAT models is non self averaging, excluding possibly the random bond 
Ising model, possibly, the theory needs some changes in order to be applied to the specific 
heat C which diverges logarithmically (and as a double logarithm for the random bond Ising 
model) and not with a simple power law. 
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^- 1 10 100 

FIG. 15. The variance of the specific heat, Vc as a function of logZ for all critical models, 
Co.. 4, and A%, B2 of the RBAT model. The solid lines are fits to the form ( |43| ) , yielding estimates 
for the fitting parameters a v ,b v , c v which are listed in table [V| . 

In fig. [16| the variance of the magnetization, Vm of the seven critical RBAT models is 
fitted by the function 

V M = a v r 2 » ln[l + c (Z {c ^ )puro - 1)] + b v l~^{\n[l + c (/ (a/!y)purc - 1)]}5 + 

c v r 2 »-i{\n[l + c (l {aMp - - 1)]}* . (44) 

The fitting coefficients a v ,b v and c v are given in Table [V]. The data are much more noisy 
than the data of the susceptibility (see Fig. |8|). 
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FIG. 16. The variance of the magnetization, Vm as a function of log / for all critical models, 
Co.. 4, and A%, B2 of the RBAT model. The solid lines are fits to the form ( |44| ) , yielding estimates 
for the fitting parameters a v ,b v , c v which are listed in table [V| . 

To summarize, we have examined the sample to sample fluctuations in various thermo- 
dynamic quantities of some random bond Ashkin Teller models. These include the random 
bond Ising and four state Potts models. It was found that far from criticality all thermody- 
namic quantities examined are strongly self averaging (that is their variance scales as l~ d ) . 
At the critical point we found that the susceptibility \i the susceptibility of the polarization 
and the magnetization M are non self averaging, while the energy E is weakly self 
averaging. The data for the variance of the specific heat seems to imply weak self averaging 
of the specific heat. Since the data are not accurate at the larger sizes used, this may well 
be a transient behavior, compatible with our theory which predicts that asymptotically the 
specific heat is non self averaging. A phenomenological finite size scaling theory was devel- 
oped for the sample to sample fluctuations. Its main prediction is that when the specific 
heat exponent a < (a of the disordered model) then, for a quantity P which scales as l p at 
criticality, its variance Vp will scale asymptotically as l 2p+ ~ . The theory is not applicable in 
the asymptotic limit (/ — > 00) to cases where - = + . Nonetheless in the accessible range of 
lattice sizes we found very good agreement between the theory and the data for V x , V X ( P ) and 
Ve- The data for V x is especially convincing. The theory also describes well the variance of 
models with weak disorder, exhibiting slow crossover to the randomness dominated behav- 
ior. The theory may also be compatible with the data for Vm and Vc, but evidence for this 
is less convincing. We note that if our assumption ( pop is incorrect and should be replaced 
asymptotically by 5T C {1) ~ l~ yt , then our theory predicts that Vp ~ l 2p independent of 
a. In this case all quantities (excluding the energy which has a non vanishing non singular 
part) are non self averaging independent of a. 

In order to further test our theory we intend to study the sample to sample fluctuations in 
the site dilute three dimensional Ising model where - < and our analysis holds. 
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APPENDIX: 

In this section we draw some more conclusions based on the finite size scaling form 
without making any assumptions on the explicit x dependence of the coefficients of Q x . 

What can we deduce about the coefficients of Q x from the Brout argument? Consider 
(p4|) in the limit I — > oo, t x finite. In this case £ is finite, the Brout argument holds, and 
one expects P X {T) to be sample independent. This means that in this limit we expect the 
coefficients of Q x and its argument to converge to some x independent values. It follows that 
we can assume that these coefficients are distributed according to some unknown distribution 
function whose width w(l) depends on I and tends to zero as / — ► oo. 

Is there any limit in which one may recover the usual finite size scaling behaviour, 
completely independent of the specific sample x7 Consider the limit I large but finite and 
T m T c . Let us now add the assumption that the width w(l) tends to zero no slower than 
l~ d / 2 . Then for laree enough I, according to equations (|20| - p3|) , t x approaches i given by 

• T ~ T c il) 

t = T ( A1 ) 
J- c 

as Z~ d/2 , and 

Q x (U Vt ) - Q(il yt ) , (A2) 



so that we recover the usual |JT| finite size scaling behaviour. 

The limit ( |A2[ ) cannot account for the large sample to sample fluctuations that we have 
numerically observed at T = T c even for rather large values of I. Indeed, special care is needed 
in the case T = T c , where £/l is not small and the Brout argument does not hold. It turns 
out that in this case the limit, where the x dependence of t x can be ignored as in ( |A"2j) , does 
not occur or is reached 'slowly'. When T = T c , then t x = —t c (x, I) = {t c (l) — t c (x, I)} —t c (l), 
so that according to (p0|) and (^) t x is a difference of a fluctuating term of order l~ d l 2 and 
a term of a constant sign of order l~ yt . Therefore 

\5i x \ f pt- d / 2 = \% = J-I&l if d/2 >y t {a< 0) 



i x \ I 1 if d/2 <y t (a> 0) 



(A3) 

Thus for T = T c and positive a, t x -/->■ t for large I and (|A2|) is not justified. In this case 
the relative fluctuations in the argument of Q x are of order 1 and their absolute magnitude 
scales as so that it increases with I. So for large I the argument of Q x is a constant plus 
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a large fluctuating quantity which increases with I. It follows that Q x cannot be expanded 
as is done in Sec. [V], and that the limit ( |A2|) does not exist. 

For negative a it follows that the fluctuations in the argument of Q x scale as Z~'^L Since 
we have assumed that the fluctuations in the coefficients of Q x , w(l), scale as l~ d ^ 2 then if 
\£\ < | then at T = T c w(l) decreases faster than the fluctuations in the argument of Q x . 
Thus one may consider the range of I for which 

Q x (U yt ) -> Q(tJ vt ) , (A4) 

where only the argument of Q x is x dependent. In this case the coefficients of Q are some 
constants for which we need not assume anything about their x or I dependence. Consider- 
ation of the limit ( |A4j) suffices to reach our main result (|T8|) (but not corrections to (jTi 



independent of the assumptions made in Sec. on the x dependence of the coefficients of 

Qx- 
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